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Abstract: Recent work in inter-jet energy flow has identified a class of leading logarithms 
previously not considered in the literature. These so-called non-global logarithms have 
been shown to have significant numerical impact on gaps-between-jets calculations at the 
energies of current particle colliders. Here we calculate, at fixed order and to all orders, 
the effect of applying clustering to the gluonic final state responsible for these logarithms 
for a trivial colour flow 2 jet system. Such a clustering algorithm has already been used for 
experimental measurements at HERA. We find that the impact of the non-global logarithms 
is reduced, but not removed, when clustering is demanded, a result which is of considerable 
interest for energy flow observable calculations. 
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1. Introduction 


The identification of the inter-jet energy flow as an infrared safe way to study gaps-between- 
jets processes has produced extensive interest in the last few years. By considering such 
observables we may start to formulate a perturbative approach to the description of such 
cross sections, as well as probe the interface between perturbative and non-perturbative 
physics at energy scales ~ IGeV. The inter-jet energy distribution was recently calcu¬ 
lated by Sterman et al ^ Q by separating out the primary emission Bremsstrahlung 
component and calculating this quantity to all-orders for a 4 jet system. However it was 
recently pointed out by Dasgupta and Salam i§ that this procedure does not include 
the effects of so-called non-global logarithms, a set of leading logarithms which were shown 
to be numerically important at the energy scales probed by current colliders. 

Recently the HI and ZEUS collaborations ii performed improved gaps-between-jets 
analyses in which the entire event is clustered into (possibly soft) jets using the inclusive 


kt algorithm [12, 13, 141. A gap event is then defined by the total minijet energy in the 


inter-jet region, rather than the total hadronic energy as in previous analyses. The hope 
was that this would ‘clean up’ the edges of the gap and make this observable less sensitive 
to hadronic uncertainties. 

In this paper we show that by demanding the gluonic final state to survive a clustering 
criterion the effect of these logarithms is reduced, but they are still numerically important, 
at HERA and the Tevatron. The study presented here is particularly interesting in the 
light of the recent analyses by the HI and ZEUS collaborations. 

The Bremsstrahlung component of the energy flow observable has been the focus of 
intense work in the last decade, for examples see i i 0, 0- It was found that by 
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considering the emission of soft, wide angle gluons by light-like quarks one may describe 
the soft gluon dynamics of hadronic processes by an effective theory known as eikonal 
theory. This subject is now very well developed and the resummation of such primary 
logarithms has recently made contact with high-but-fixed order calculations |^]. This 
formalism now allows detailed calculations of inter-jet energy flow as well as illuminating 
insights into the topology of colour mixing ||^. 

However, in the study of single-hemisphere observables |^] Dasgupta and Salam iden¬ 
tified that an important class of logarithms was missing from the soft gluon calculations. 
These contributions arise in observables that are sensitive to radiation in a restricted re¬ 
gion of phase space and are formally the same order as the primary emission component 
for inter-jet observables. To see the origin of these logarithms more clearly, consider soft 
gluonic radiation into a patch of phase space H arising from a 2 jet system, where we re¬ 
strict the total energy of radiation into H to be less than Qo. The primary logarithms arise 
from gluons that are emitted directly into H with energy vetoed down to a scale Qn and 
form the single logarithmic (SL) set (a^ log ), n> 1, where Q is the scale of the jet 

line. Now consider a gluon being emitted outside of H with intermediate energy Qi, and 
then vetoing subsequent emission from this gluon into H down to scale Qq. Integrating Qi 
up to Q then generates another SL set of (a^ log ), n > 2, which are formally the 

same order as the primary emission terms. This was studied in detail in [^, and in this 
work we consider how the conclusions found are modified phenomenologically when the kt 
clustering algorithm is applied to the final state. 

The organisation of this paper is as follows. Section 2 describes the kt clustering algo¬ 
rithm used in the HI and ZEUS analyses, known as the inclusive jet clustering algorithm. 
We then develop the algorithm, in order to derive a form suitable for use with a 2-gluon 
system. Section 3 then describes the order calculation of the effect of non-global loga¬ 
rithms in our 2-jet system, with the clustering algorithm included in the calculation. We 
find that the asymptotic suppression factor with clustering is reduced with respect to the 
non-clustering case. We then proceed to calculate the non-global effect to all orders in 
section 4 using the combination of the large W limit and a Monte Carlo algorithm and we 
conclude with a summary and a discussion of further work. 

The all-orders treatment used in this work allows direct inclusion of the kt algorithm, 
in exactly the same way as is used experimentally, and we find that the clustering process 
reduces the magnitude of the non-global corrections to the primary suppression factor. 
However they are still phenomenologically relevant at HERA. 


2. The kt clustering procedure 

The version of the algorithm we use, which is the one used in the HERA analyses, is 


known as the inclusive kt algorithm [12, 13, 14|. The main features of importance to 


the present study are: the clustering procedure starts from the particles of lowest relative 
transverse momenta and iteratively merges them to construct pseudoparticles of higher 
transverse momentum; the decision of whether a particular pair of pseudoparticles are 
merged depends only on their relative opening angle (see below); despite this, it is possible 
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for soft particles to be ‘dragged’ through relatively large angles by being merged with harder 
particles, which are merged with even harder particles, and so on. Nevertheless, we do not 
expect out results to be qualitatively different from other infra-red safe jet algorithms, 
such as the improved Legacy Cone algorithm of |17]. In the remainder of this section we 
will provide a formulation of the kt algorithm that can be directly applied to our fixed 
order calculation in the next section. We follow the HI and ZEUS analyses closely and 
set the radius parameter, R, to unity for most of the numerical results. The algorithm 


is implemented in the package KTCLUS []^|, which is used for the all-orders calculation 
later in this work. 

It is necessary to use the full iterative algorithm for experimental analysis and for 
Monte Carlo applications, but when we consider a 2-gluon final state in this work we can 
reduce the algorithm to a convenient analytic form. We start by considering a hard jet 
line at scale Q which radiates a gluon with some transverse energy Et^i in some direction 
(7/i,(/>i). This system then radiates a secondary soft gluon with transverse energy Et ,2 in 
some direction ( 772 , 4>2)- We assume that the transverse energies of the gluons are strongly 
ordered, 

Et,1 Ej'^2- (2T) 

The kt algorithm works by defining for every pair of particles ij a ‘closeness’ dij and for 
every particle i a ‘closeness to the beam direction’ di. For the 2 gluon final state, these are 

di = E‘^1^ 

d2 = 

712 


di2 — Erp2[{rji — 772 ) -|- {(l>i — ^ 2 ) ]■ 


( 2 . 2 ) 


By considering the strong ordering of the transverse momenta, the two gluons will be 
clustered if dij < d 2 , so we require 


(f?i - V2f + {4>i - hf > R? 


(2.3) 


for the two gluons to constitute separate jets and not be merged by the algorithm. Therefore 
for a two-gluon final state to pollute the gap and generate secondary logarithms we require 
that gluon 1 is outside the gap, gluon 2 is inside the gap and that they be sufficiently 
separated in ( 77 , </>) to avoid being merged. Therefore the clustering condition manifests 
itself as a ©-function in our calculation, 

+ (2.4) 

where we have used our freedom to set = 0. We will use this result in the next section. 


3. Fixed order calculation 

3.1 Definitions and primary emission form factor 

Following Dasgupta and Salam Q, the observable we are interested in is the total transverse 
energy Et flowing into a region of phase space H for an event characterized by the hard 


- 3 - 




scale Q, 


(3.1) 


Et — Et^i- 

ien 

We are specifically interested in the cases of Q being either a slice in rapidity or a patch, 
bounded in rapidity and azimuthal angle. The quantity we shall calculate is called Sq and 
is defined to be the probability that Et is less than some energy scale Qq. 




1 

Jo 


rQn da 
dEt-r-^ 
dEf 


(3.2) 


We shall assume the strong ordering <C Q. The aim of this work is to calculate the 
importance of the non-global contribution to Sq and so it is convenient to factorize this 
expression into a function describing primary emission into fl, Sn,p(t), and a function 
describing (secondary) emission into from large-angle soft gluons outside of fl, 5(t), 


Ho(t) = S{t)Efi,p{t). 

(3.3) 

We have denoted the following integral of ag by t. 


1 /■'5/2 d kt 


t{Qn,Q) = — / ^asikt), 

(3.4) 

1 , f oisiQn) \ 

~ AnPo °"U.(Q/2)y ’ 

(3.5) 

0's , Q 

~ 27r 2Qn ’ 

(3.6) 


where the first equality is exact, the second holds at one loop, the third assumes a fixed 
coupling and /3o = (IIC^ — 2nj)/(127r). The leading order contribution to S{t) comes in 
at and we shall calculate this for a 2 jet system in the next section but it is useful to 
first consider the primary emission function at first order in a^. If we do not restrict the 
phase space for gluon emission then, order by order in perturbation theory, we expect a 
complete cancellation of real and virtual soft gluon contributions to the primary emission 
form factor. However the requirement of a gap in a restricted region of phase space results 
in this cancellation being spoiled and we are left with an integral over the vetoed region. 
Hence, to order as, 


E^n\Qn,Q) 



(3.7) 

(3.8) 


where Aq, denotes the area in (ry, 4>) space of the region H. The assumption that the primary 
gluons are emitted independently according to a two-particle antenna pattern means that 
by exponentiating the one loop answer and running the coupling to the scale kt, we can 
write T,Q^p{t) to all orders, 


^n,p{t) = exp {-ACpAnt). 


(3.9) 


This equation only includes contributions from independent primary emission. 
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3.2 The function S{t) 

The non-global contribution to Sq is contained in the function S{t), which has its first 
non-trivial term at order Hence we can write the following expansion for S{t), 

S{t) = l + S2t^ + S3t^ + ... = l + Y^ Snt^. (3.10) 

n=2 


The goal of this section is to calculate S 2 for two different geometries of H with the condition 
that the topology of the gluon tree satisfies the clustering algorithm defined previously. We 
begin by defining the following 4-momenta, 


h 

h 

k2 


|(1,0,0,-1), 


(3.11) 

|(1,0,0,1), 


(3.12) 

^t,i(cosh(r/i),sin(()>i) 

,cos((/>i),sinh(r/i)), 

(3.13) 

Pt,2(cosh(r/2),sin(())2) 

,cos(y)2),sinh(r/2)), 

(3.14) 


and for a general region H define S 2 by 

(4i) {£)) - 

[ / X 2 dx 2 [ xidxiQ fx 2 - Ws, (3.15) 

Jk2€n 27r 16 Jq \ Q J 

where we define the transverse momentum fraction by 

Pt,i=xi^. (3.16) 

This expression for S 2 contains the secondary part, Ws, of the well-known matrix element 
squared for the energy ordered emission of two gluons, which can be derived from gluon 
insertion techniques [p^], 

W = Ar (CaJ^ Ca (61) / _ {ah) \ 

^(al)(16) V 2 (a2)(21) 2 (62)(21) V ^ 2 y {a2){2h)) 

= ClWp + CfCaWs, (3.17) 


where the notation {ij) denotes the dot product of the appropriate 4-momenta. This 
expression contains the primary emission piece Wp, proportional to Cp, and the piece that 
interests us, which is the part proportional to CfCa and denoted Ws- Note that the last 
term is the dipole interference term and is absent in the large W limit. Denoting the 
secondary emission piece by W 2 and evaluating the 4-momenta products we get 

{ab) = (3.18) 

(la) = ^^exp(-77i), (3.19) 
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(3.20) 


(16) = 
(2a) = 
(26) = 
( 12 ) = 


Q'^xi 


Q‘^X2 


exp(+ryi), 


exp (- 772 ), 


Q‘^X2 
4 

Q‘^XiX2 


exp(+r? 2 ), 

(cosh(? 7 i - r/ 2 ) - cos((/>i - 4 ) 2 )). 


Now setting (j)i equal to zero we arrive at 

128 f cosh(7/i — 772 ) 




Q'^x'lxl Vcosh(77i - 772 ) - cos(())2) 


- 1 


(3.21) 

(3.22) 

(3.23) 


(3.24) 


The energy fraction integrals are straightforward and lead to the following leading logarithm 
(LL) result (the 0-function ensures that sufficient energy reaches the gap region Q,) 




Hence we arrive at 


52 = -ACfCa 


L 


angles 


cosh (771 - 772 ) 


cosh(77i — 772 ) — cos{4>) 


- 1 


(3.25) 


(3.26) 


Specializing to a slice in rapidity of width A 77 and delimited by I 77 I < A 77 / 2 , where boost 
invariance allows us to center the gap on 77 = 0 , we can exploit the symmetry of the 771 
integral and perform the trivial (p integral of the first gluon to express the final integral 
over angles as 


L 


angles 


At] 


= 2 


dr]i 




d4>2 

27r ' 


(3.27) 


For the case of no clustering of the gluons it is possible to perform the azimuthal average 
using the result, 


p27r 

Jo 


dp 


1 


1 


(3.28) 


Iq 27r cosh(77i - 772 ) - cos((/>2) | sinh(77i - 772 )] ’ 

which can be proved by noting that the following contour integration over the unit circle, 


I=- 

i 


d z 


(z - exp(7j2 - r]i)){z - exp(77i - 772 )) ’ 


(3.29) 


is equivalent to the integration over the azimuthal angle of the second gluon. For the case 
of no clustering of the gluons, it is possible to perform all of the integrations analytically 
and obtain |^, 


52 


-4CfCa 


12 


+ (At?)^ 


A77log(e^^^ 


1 ) 


- iLi 2 (l - e“”) . 

(3.30) 
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(3.31) 


lim 

Atj—^oo 


= -CfCa- 


Note that as At/ increases, S 2 rapidly saturates at its asymptotic value, 

2^ 

3 ■ 

This analytic evaluation of S 2 is not possible for the case of clustered gluons because, as 
discussed previously, the requirement that a gluonic final state is in a configuration that 
will survive a clustering algorithm can be written as a 0-function of all three angular 
integration variables, 

+ (3.32) 

However we can readily reduce the three-dimensional integral to a one-dimensional integral 
using standard techniques if we consider the region of phase space that is vetoed by the 
clustering algorithm, denoted S^- Doing this we obtain, 

-32CfCa 


S^2 = 


TT 


f 

/ min(r/, Arj) 

Jo 


2 coth(r/) arctan 


( tan(v^i?2 _^2/2) \ ^ ^ 

tanh(T//2) J ^ 




dr], 
(3.33) 

where we define 7 / = t/i — 7 / 2 . Therefore the solution for S 2 with the clustering condition 
imposed is obtained by subtracting from the analytic unclustered result. We resort 
to numerical techniques to solve this equation, which has the advantage of being easily 
extendible to any gap geometry of interest. Our numerics were done using Monte Carlo 
integration methods and figure shows our result for S 2 as a function of At/ for a cluster 
radius oi R = 1.0. In this figure we also show, as the solid line, the curve for S 2 which is 
obtained with no clustering condition imposed on the gluons. In |p the saturation for S 2 
was explained by observing that the dominant contribution to S 2 arises from the situation 
with one gluon just outside the gap emitting a gluon just inside the gap. This occurs when 
7/1 ~ 7/2 — — A 7//2 and the dominant contribution arises from the collinear region of the 
matrix element. The results we have obtained show that when we demand the gluonic final 
state to survive a clustering algorithm, the saturation of S 2 observed in Q is still observed 
but the saturation value is reduced by 70%. In other words the value that S 2 saturates to 
is reduced from 6.57 to 1.81. The reason for this reduction is that we have removed gluons 
from the region of collinear enhancement but gluons that are still sufficiently separated in 
the (7/, (/)) plane to satisfy 

{m - + 4'2> (3-34) 

will survive the clustering process and contribute to the non-global logarithms. Therefore 
we conclude that the saturation of the non-global contribution at fixed order is still seen 
when we demand clustering of the final state. If this reduction holds to all orders then the 
magnitude of the non-global logarithms will be smaller for the clustered case relative to 
the non-clustered case. To explore this it is necessary to perform the all-orders treatment, 
which we do in the next section. 

4. All-orders calculation 


The extension of the above calculation to all orders presents considerable mathematical 
problems due to the geometry and colour structure of the gluon ensemble. Therefore we 
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have used the large approximation and numerical methods, as developed by Dasgupta 
and Salam |^, to extend the calculation to all orders. We have extended the simulation 
developed in the work of Dasgupta and Salam to include clustering of gluons according to 
the kt algorithm discussed previously. This is performed by terminating the gluon cascade 
evolution if D is polluted and the final state survives the clustering algorithm. If the latter 
condition is violated, the Monte Carlo algorithm is allowed to continue. 

Note that this algorithm can produce configurations with many gluons at high t and 
this, coupled with the fact that the speed of the kt algorithm scales like (where N is 
the number of gluons), means that the Monte Carlo algorithm can be very slow in this 
region. The result is that we have large statistical errors for t > 0.3, so we do not show 
this region. In the next section we shall present our results. 


4.1 All-orders results 


We have verified the results obtained by Dasgupta and Salam Q, which have been calcu¬ 
lated with no clustering requirement imposed on the gluons. 

Figure § shows the function S{t) for two different geometries for Q: a slice in rapidity 
of width Ar] = 1.0 and a square patch in rapidity and azimuthal angle of side length 
Ar] = A(p = 2.0, with the requirement of clustering on the final state gluons. This is 
implemented in the Monte Carlo algorithm in the same way as it is used by HI and ZEUS, 
allowing us to estimate the experimental impact of the clustering scheme on the non-global 
logarithms. Firstly the figure shows that at high t {t > 0.2) the suppression increases 
faster for the full calculation than for the exponentiation of S 2 with clustering and the 
two curves have different shapes in this region. This agrees with the analysis of the case 
with no clustering. Secondly we see that the t dependence of the suppression is geometry 
independent at high t; the implication of this is that the clustering process maintains a 
buffer region^ of suppressed intermediate radiation around D exactly as in the unclustered 
case. Therefore figure |2| shows that the clustering requirement does not change the gross 
features of the non-global suppression at all orders and the results indicate that the buffer 
region postulated to exist around the gap D is preserved by the clustering process. 

Figure ^ shows the reduction of the phenomenological significance of the non-global 
logarithms when we cluster the final state. The figure shows the function S(t) with only 
primary emissions and the full all-orders treatment with and without clustering effects. 
This is done with D defined as a slice in rapidity of size Arj = 1.0. There are several 
points to note. Firstly the effect of the non-global logarithms is a large suppression of the 
cross section relative to the primary-only result. The value of t which we can consider the 
realistic upper limit for next generation experiments is about 0.15 so we shall take t = 0.15 
as our reference value. We can translate this into a energy scale by using the running 
coupling definition of t, equation and obtain Q ~ 100 GeU for Qq = IGeV. In this 
region the inclusion of the non-global effects without clustering increases the suppression 
relative to the primary-only result by about a factor of 1.65. When we include clustering 


^The buffer region is defined by the absence of any reconstructed jets within it. It may therefore contain 
gluons, provided that they get pulled out of the buffer region by the clustering algorithm. 
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effects, this difference is reduced to about 1.2. Therefore, at all orders, the requirement 
of clustering on the final state reduces the phenomenological significance of the non-global 
effects by about 70%. This reduction in magnitude of the effect can be seen to persist 
to all orders. Hence when calculating cross sections, if we exclude the effect of non-global 
logarithms then we will overestimate the cross section by 65% for a non-clustered final state 
and by 20% for a clustered final state. At larger t values, the overestimation increases. For 
comparison, the typical errors on the HI gaps-between-jets data is ~ 30%. 

Finally, figure ^ shows the same comparison of the full result with and without clus¬ 
tering for a patch of size Arj = A(j) = 2.0. We can see that the conclusions we made for 
figure ^ apply and the effect of the non-global logarithms is of similar magnitude. 

We have also performed calculations to see how the non-global suppression is affected 
by varying the radius parameter, R, of the clustering procedure. By decreasing R the 
impact of the clustering algorithm is reduced and the effect of the non-global logarithms is 
restored to the non-clustered case. We expect the magnitude of the non-global suppression 
to tend to the non-clustered case as 7? —> 0. Similarly, increasing the radius causes more 
gluons to be included in the clustering and hence the magnitude of the suppression to 
reduce. In fact, we have found that for R close to 1.5 the full result is almost identical to 
the primary result. 

In summary, non-global logarithms are important not only from the point of view of 
correctness of the leading logarithm series but also result in significant numerical corrections 
to cross sections. These corrections are reduced by about a factor of 3 if we cluster the final 
state. It is clear, therefore, that while the use of the kt algorithm has aided the control of 
the non-global logarithms, they still have a significant numerical effect at HERA. 

5. Conclusions 

The observation that observables that are sensitive to radiation in only a restricted part of 
phase space, so-called non-global observables, are strongly sensitive to secondary radiation 
is a new and exciting discovery. For a long time it was widely thought, now it seems 
incorrectly, that it was sufficient to only consider primary emission contributions to such 
observables. These primary, or Bremsstrahlung, contributions are well understood for a 
variety of processes. Recent measurements of gaps-between jets events at HI and ZEUS are 
exactly the class of observable that are sensitive to these effects and to deal with this fact, 
a study of non-global logarithms in the context of HERA measurements is required. In 
this work we consider the non-global logarithms generated from secondary radiation into a 
restricted region of phase space under the condition of the final state surviving a clustering 
algorithm. Such kt clustering algorithms have been used for the HERA measurements. 
In this paper we have extended the study of Dasgupta and Salam and found, for a 2-jet 
system with clustering imposed, that the clustering process reduced, but did not eliminate, 
the non-global logarithm effect. 

Our main conclusion is that the final state specified in the HI and ZEUS analysis will 
mean that a primary emission calculation in the manner of Sterman et al i i 0 I of 
energy flow observables will overestimate the observed gaps-between-jets rate by around 
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20%. This is to be contrasted with an overestimation of 65% that would be found for a 
non-clustered final state. This result was calculated to all orders in the large Nc limit for 
a 2 jet system. We reserve the extension of these results to the more complex 4 jet case 
for future work but we expect the gross conclusions of this study to apply. 
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Figure 1: S '2 as a function of Ary. The clustered curve is obtained using a kt algorithm with a 
radius parameter of 1 . 0 . 



Figure 2: The function S{t) for a slice and a patch of phase space, with the condition of kt 
clustering imposed on the gluons. These curves were obtained using a Monte Carlo procedure in 
the large limit. Also shown are the curves for S(t) obtained by exponentiating the fixed order 
result. The geometry independence of the t dependence indicates that the buffer zone mechanism 
identified to exist in previous work survives the clustering algorithm. 
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Figure 3: The phenomenological implications of clustering for the case of ft being a slice of rapidity 
of width 1.0. The reduction of the non-global suppression factor when clustering is included can be 
seen. 



Figure 4: The effect of the clustering of the gluons on the function S{t) for a patch of phase space 
of size Arj = Acj) = 2.0 
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